A MATHEMATICAL APPROACH TO THE 
PLATO'S PROBLEM 
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ABSTRACT. Maybe the first inverse problem presented in the history of the occi- 



dental thought is described in the book Republic, written by Plato. The problem is posed 
in the Book VII in a text known as the Allegory of the Cave. That text motivated us to 
formulate a simple mathematical model that simulates, in a sense, the situation of the 
persons described in that problem. 
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\Q , 1. INTRODUCTION. The motivation of this paper is a problem posed more 

Os 

than 2000 years ago by Plato. Nevertheless the essence of that problem is one of the most 

■ challenging questions to be found in modern science and technology. Plato in the Book 

O. 

VII of his monumental masterpiece Republic, wrote a dialogue known as the Allegory of 

O ■ 

the Cave [lj . That text is maybe the first inverse problem presented in the history of the 

•rH , 

^ ■ occidental thought. Briefly, several persons, sitting, facing the rear wall of a cave were 

totally immobilized even unable to move their heads. A fire was set between the persons 
and the wall. A procession of men and women of all ages and animals carrying, pulling 
and pushing all kind of objects moved between the wall and the fire. The persons tied 
on the chairs were able only to watch the images projected on the wall but could not see 
the parade directly. The challenge was to interpret the images and find out what kind of 
persons, animals and objects composed the parade. Plato explored through this allegory 
the knowing process and how to reach the truth. Despite the fact that his aim was distinct 
he presented maybe for the first time in a well organized structure an inverse problem, 
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specifically, pattern recognition. The Allegory of the Cave is one of the most marvelous 
texts ever produced by the human mind. 

Plato's text motivated the design of the kit presented in the sequel that simulates in 
a simplified version the images observed by the persons tied to the chairs in Plato's story. 
The first part of this paper will describe the direct problem and the last section will deal 
with a particular inverse problem. 

Nowadays one of the most important challenges that falls into this category is the 
pattern recognition problem. There are several applications in medical diagnosis through 
the interpretation of images obtained with sophisticated equipments as computerized 
tomography and nuclear magnetic resonance devices. The question is: given the pattern 
what is the tissue density and distribution in the organ under examination. 

The problem posed here is essentially the same except that the number of parameters 
involved is substantially smaller. If we know how to solve the direct problem the solution 
of the corresponding inverse problem can be substantially facilitated. For the reader 
interested in more information on inverse problems see [2]. 

2. THE PROBLEM. A slender object with the shape of a thin plate B rotates 
with angular velocity Q around an axis O perpendicular to a plane containing the object. 
Such an object is confined inside a box C such that the slot F and the object belong to 
the same plane. 

A device inside the box, consists of a lamp and a concave mirror which projects the 
profile of the plate B on a moving screen E. This projection produces a shadow with 
height H equal to the distance between the upper and lower tangents to the contour of 
the plate B. 

If the film is sensitive to light, and unwinds in front of the slot with velocity v at the 
same time that the plate rotates around O, an image which is in a certain way a "picture 
of B" will be generated. 
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Figure 1: Building setup: object and projection system 
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Figure 2: Image of the object B recorded on the film 

The image printed on the tape is a continuous picture with the shape of a ribbon 
with variable width, we will call it "kinematic image" of B. 

The picture is bounded above and below by Y a (z) and Y^z), respectively. The width 
is given by H(z) = Y s (z) - Y^z). 

The following problems arise: 
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i) Given the geometry of B, the angular velocity Q(t) and the travelling velocity v(t) 
find H(z,to) = Hq(z) after a given time to- 

ii) Given Y s (z),Yi(z) and H (z) find the geometry of B. Is the solution unique? 

The first problem is known as a "direct problem" and the second one is an "inverse 
problem" . 

Observations: 

1. The projection Hq(z) do not identify reentrancies. In other words, it is a critical 
condition the contour to be convex in order to have an unique solution. 



2. Another difficulty is the existence of vertices on the contour. These singularities 
could be overcome, but they will need special treatment since at the vertices there 
is not a definite tangent. We will see that for the case of regular polygons it is 
possible to determine a simple solution. 




Figure 3: Non-convex contours 




Figure 4: Contour with vertices: tangents not well-defined 
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3. FORMULATION OF THE DIRECT PROBLEM. Consider a flat object B 
with a smooth contour. Let the contour, be defined in terms of polar coordinates r(/3)) 
having as a reference Qxy fixed in B. The object B moves with respect to OXY rotating 
around a perpendicular axis to its plane, going through the point O fixed in OXY, with 
angular velocity 9'{t) = d9/dt, where 6{t) is the angle between a line Qx fixed in B and 
a line OX fixed at the inertial reference system. 
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Figure 5: Contour and tangent representation 



Our propose is to find the shadow of B projected on the rotating tape, that is, the 
distance H(z, t) between the upper and lower tangents. 

Let f(P) be the position vector of an arbitrary point M at the contour of B. We 
are admiting a smooth contour, that is, without corners. The general expression of the 
tangent related to the contour of B is: 

df 
1 = df3~ 
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Figure 6: Limits Y s and Yj of the projection and height H 

Clearly the upper and lower tangents are paralell to the horizontal axis OX. This is 
equivalent to say that both are orthogonal do OY. Therefore we may write: 



t-j = o, 



(3-1) 



where j is the unitary vector corresponding to the axis OY. 

The representation of t in OXY is given by T = IZt where the rotation matrix is: 



K 



cos 9 — sin 9 
sin 9 cos 9 



The condition (13. ip can therefore be written as: 



Kb ■ j = 0. 



(3.2) 
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If the contour is smooth and convex there will be two solutions for this equation 
corresponding to the upper and lower tangents. This expression relates (3 and 6 in such 
a way that with this relationship (3 will always define the same position of the extreme 
tangents. We will designate this position (3(6) = (3* . 

In general there is not an explicit form (3(6), and hence, (3* has to be determined 
from the numerical solutions of (13.21) . In the next sections we will to study some special 
cases where it is possible to obtain (3* explicitly. 

Once we determine (3*, or (3* and (3* corresponding to the upper and lower tangents 
it is possible to find Y s and Yf. 

Y s = R(P* u ,6)-j = R u -l 
Y l = R(P*,6)-j=R r j, 

where R u and Ri are position vectors defining the tangent points M and N on the 
contour corresponding to the upper and lower tangents represented in OXY, that is: 

R u = K(r Q + r((3*)), 
R l = Tl(f Q + ?(/%)). 

Now, 

H = Y S — Yi = [n(f((3* s ) - r(/3* s ))} ■ J, 

that is, the total height of the shadow does not depend upon the pole of rotation O. 
Example: The figure 7 shows two projections of the same circle with different poles. 

It is now necessary to find the image printed on the moving film that travels in front 
of the small hole with velocity v(t). 

Consider the axis O'Z on the film, orthogonal to the OXY plane. Let O' be the 
origin of the Z axis. Clearly the figure printed on the film, that is, the kinematic image 
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Figure 7: Image of a circle. Rotation around the center and around a point on the contour 

of B, will be comprised between Y s and Y$, corresponding to the points where the upper 
and lower tangents intercept the tape. 

The projection Y s on the film varies with time as: 



(TY d f) df) 



On the other hand the film moves with speed v, therefore 

dz , . 

* = "<«>■ 



The solution of the system 



dY s 
~W 

dz 
dt 



v(t) 



gives Y s = Y 8 (z). Similarly we would obtain Y { = Yi(z). 



Note that the solution of the equation (13.21) has two non-coincident roots if and only 
if the object B is convex and its contour is a smooth curve. 
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Figure 8: Non-convex contour with more than two tangents 



4. THE CASE OF THE CIRCLE. Let a circle of radius a rotates around the 
centre O. 
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Figure 9: Image of a circle rotating around the center 



The solution is trivial. The kinematic image is a strip with width equal 2a, bounded 
by Y s = a and Y{ = —a. Formally: 





cos (5 




— sin P 


r = a 




,t = 






sin p 




cos P 



Kt = 



cos 9 — sin 9 
sin 9 cos 9 



— sin /3 

cos/3 



- sin(fl + /3) 
cos(# + P) 
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The condition ( 13.2ft TZt ■ j = implies that 008(6* + /3) = 0, that is, 



7T 



+ /? = (2n-l)-, n = l,2. 



So, we have two different solutions: 

fts = \ - 9, 
Pi = tj- — 0. 

Suppose that the disc with radius a rotates with constant angular velocity 9 
fore it is possible to write: 



ut. There- 



3tt 



cut. 



Now, the general expression of the position vector on OXY is 



R = a 





cos ft 




cos{9 + ft) 






= a 






sin ft 




sin(fl + ft) 



It follows that 



cos 9 — sin 9 
sin 6* cos 9 



Y s = R(ft s )-j = asm(9 + ft s 



Since /3 S 



^ — 9 it comes: 



Y» = a sin 



7T 



Similarly, Yi = —a and = Y s — Y> — 2a. 



A more interesting case arises when the disc rotates around a point that does not 
coincides with the center. Take the case of a disc rotating around a fixed point on the 
contour: 



r 



1 




, 9(t) = cut. 
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Figure 10: A circle rotating around a point of the contour 



There is no loss of generality if we take the x axis fixed on the disc through the 
rotation center. Following the reasoning exposed above, we have successively: 





cos (3 




— sin P 


r = a 




,t = 






sin p 




cos P 



nt 



- sin(fl + P) 
cos(6 + p) 



The condition f)3.2p IZt ■ j = implies that 



Ps 

Pi 



= 5- 0, 

3n a 
-T~ 9 
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The general expression of the position vector on OXY is 



R = TZ(f + r) 



cos — sin 9 
sin 9 cos 9 





( 


i 




cos f3 






a 




+ a 




) 




K 







sin j3 





Its projection on the Y axis is 

R-J=a [sin 9 + sin(# + (3)] = Y. 

Substituting (3 by (3 S and $ we find 

y s = y(/3 a ) = a[sin(0) + l], 
= = a[sin(^) - 1] 

and // = Y s - Y { = 2a. 

The printed image over the film that winds down in front of the slot with velocity 
v = v o =constant is obtained from the system 



^= acos^f = acucos(cjt), 
dz _ ,, 

it- v °- 

Using the conditions 5^(0) = a and z(0) = the solutions of the equations above are: 

Y s = a (1 + sin(c<jt)) , 

z = Vot. 

Solving this system, by elimination of the variable t we obtain: 

Y = a ( 1 + sin ( — z ) ) . 



Similarly we also obtain: 



Yi = Y s - H = a [ -1 + sin ( — z 
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Figure 11: Kinematic image of a circle rotating around a point on the contour 
The kinematic image has the shape of a wave (see figure 11). 

5. THE CASE OF THE ELIPSE. Consider the position vector of an arbitrary 
point on the contour of an elipse centered in the origin of a system OXY with semi-axis 
a > b: 



a cos (3 



r = 



b sin (3 



Its tangent vector is given by: 



t = 



—a sin (3 



b cos/? 



The rotation of tangent vector by an angle 6 is: 



Kt 



—a cos 6 sin (3 — b sin 6 cos (3 



—a sin 6 sin (3 + b cos 9 cos (3 



The tangency condition (13.21) leads to 



a sin 6 sin (3 = b cos 6 cos (3. 



(5.1) 
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The general expression of the position vector referred to OXY is 



R 



cos 9 — sin 6 
sin 9 cos 9 



a cos (3 
h sin (3 



a cos # cos (3 — b sin 9 sin /? 
a sin 9 cos /3 + 6 cos 9 sin /? 



and its projection on the axis OY is 



Y = R ■ j = a sin 9 cos (3 + b cos # sin (3. 



(5.2) 



Rather than try to write f3 as a function of # let explore useful relationships between 
(15.11) and (15. 2p . In fact, squaring (15. 2p we have 

F 2 = a 2 sin 2 6* cos 2 /3 + b 2 cos 2 6* sin 2 /3 + 2ab sin 6* cos (3 cos 6* sin /3. 



Using the tangency condition (15 .ip we obtain successively 

y s 2 . = a 2 sin 2 9 cos 2 /? + b 2 cos 2 6» sin 2 (3 + 2a 2 sin 2 6» sin 2 (3 



,2 sin 2 + b 2 cos 2 6 1 sin 2 (3 + a 2 sin 2 # sin 2 /3, 



that is, 



Similarly we get: 



Y 2 { = a 2 sin 2 9 + (a 2 sin 2 9 + b 2 cos 2 0) sin 2 /3. 



y s 2 . = fo 2 cos 2 + (a 2 sin 2 + b 2 cos 2 0) cos 2 /?. 
Adding (15.31) and (15. 4p we conclude that 



(5-3) 



(5.4) 



which means that 



Y 2 { = a 2 sin 2 9 + b 2 cos 2 9, 



Y s = Va 2 sin 2 9 + b 2 cos 2 9, 
Yi = -Va 2 sm 2 9 + b 2 cos 2 9. 



(5.5) 



Since 9 = ut and z = Vot, the printed image on the film is obtained from the equations: 



CO 



Y s = J a 2 sin ( — z + b 2 cos 2 — z , . 



CO 



(5.6) 



Y 



a 2 sin 2 ( — z ) + b 2 cos 2 ( — z 



UJ 
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The figure 12 shows the printed image on a moving film of an elipse rotating around its 
center with constant velocity v : 

Y 




Z 



Figure 12: Kinematic image of an elipse rotating around its center 

7. THE CASE OF REGULAR POLYGONS. The techniques presented in 
the section 3 cannot be applied to this case since we are dealing with plane figures with 
corners, non-smooth contour. We will see that for this case the profile projected on the 
moving screen is defined by the vertices or by the corresponding position vectors. 

For sake of simplicity let us call P the vector OP where O is the origin of the 
coordinate system and P a general vertex of the polygon. Now given a vector P = [x, y] T 
, the Y-component after a rotation 9 is given by: 

fp(0) = x sin 9 + y cos 9. 

Let us call this function RotProj. Consider the polygon inscribed in a circle whose 
center coincides with the origin of the X — Y axis system. Assume that the polygon 
rotates around the point coinciding with the center of the circle with constant angular 
velocity u> and the film moves with constant velocity v . 

Case 1: the square. Let us consider a square with side equals a as shown in the 
figure 13. 
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Figure 13: Image of a square 
For 8 = the vertices are given by: 

A =|(1,1), 
£ = §(-1,1), 
C = f(-1,-1), 
D = f(l,-l). 

The corresponding functions RotProj read: 

§(sin(9 + cos 8), 
f (sin# - cos 8), 
|(- sind - cos#), 
§(-sin(9 + cos(9). 

Now when the square rotates of an angle 8 within the interval < 8 < n/2, the 
vertices defining the upper and bottom contours of the profile projected on the film are A 



file) = 
fsifi) = 
/<?(*) = 
hi?) = 
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and C respectively. The vertices B and D in that case are "hidden" (see figure 14). The 




Figure 14: Successive positions of a rotating square 



expressions for the upper and bottom contours are given by Y s = f^{9) = a/2 (sin 9 + cos 9) 
and Yi = —Y s respectively. 

For the case where 9 belongs to the interval 7r/2 < 9 < tt, the control vertices are D 
and B while A and C remain hidden. 

Following this procedure till the rotation completes a complete cycle (2ir) the expres- 
sion from the upper and bottom contours are given by: 



' §(sin# + cos#), if < 6» < tt/2, 

§(sin#-cos#), if tt/2 < 6> < tt, 

f(-sin0-cos0), if tt < 9 < 3tt/2, 

§(-sin# + cosfl), if 3tt/2 < 9 < 2vr, 



and Yi = —Y s . 



The image projected on the film, that is, the kinematic image, is depicted in the 
figure 15. 
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Figure 15: Kinematic image of a square rotating around its center with constant angular 
velocity. Film moving with constant speed v — v 

Case 2: The Equilateral Triangle. Consider now an equilateral triangle with side 
equal a as shown in the figure 16. 

4 Y 




Figure 16: Image of an equilateral triangle 



For 9 = the vertex coordinates are given by: 



A=^I(>/3,1), B = a^(-V3,l), C= ^(0,-1). 
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The corresponding projections are given by the functions: 

ay/3 



'(V3sin# + cos#), 



6 

f s (0) = ^(-^sin^ + cos^), 
f d (9) = -^ cos 9. 



When the triangle rotates of an angle 9 comprised within the interval < 9 < n/3, 
the upper contour is determined by A and the bottom contour by C . In that case the 
vertex B is hidden (see figure 17). 
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Figure 17: Sucessive images of an equilateral triangle rotanting around its center 



Therefore, within this angular interval the upper and bottom projections are given 

ay/3 ay/3 
by Y s = fx(6) = — —(a/3 sin ^ + cos#) and Yi = fg{0) = — cos9, respectively 

O O 

For the next interval n/3 < 9 < 2n/3 the projection is controlled by A and B with 
C hidden. 

The table below shows the intervals composing a complete cycle (2tt), the correspond- 
ing control vertices and the hidden vertex as well. 
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Angular interval 


Upper control vertex 


Bottom control vertex 


Hidden vertex 


U \ V \ 71/ O 






i > 




A 


J J 


a 


27r/3 < < 7T 


C 


B 


A 


7T < 9 < 4tt/3 


C 


A 


B 


4tt/3 < 9 < 5tt/3 




A 


C 


5tt/3 < 9 < 2tt 


5 


C 


A 



Using the table above the respective functions RotProj give the upper and the bottom 
contour curves. They are: 



a ^(v / 3sin^ + cos^), if < 9 < 2tt/3, 
6 

" V cos#, if 2tt/3 < 6 < 4tt/3, 



and 



3 

a\/3 
~6~ { 



-VSsin9 + cos9), if 4tt/3 < 9 < 2vr, 



t\/3 



3 

cza/3 



cos 9, 



H0<9 < tt/3, 



6 

a\/3 



6 

a\/3 



(-v^sintf + costf), if tt/3 < 6* < tt, 
(VSsm9 + cos9), if tt < 9 < 5tt/3, 



cos 6 1 , 



if 5tt/3 <9<2tx. 



The image printed on the film is given by figure 18. 

Exercice: Find the shape of the kinematic image on the film of a regular polygon 
having n sides. What is the difference between the image of the polygons consisting of 
even and odd sizes? 
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Figure 18: Kinematic image of an equilateral triangle rotanting around its center 



7. THE INVERSE PROBLEM. We have discussed the solution of the direct 
problem. That is, find the image of a bi-dimensional rotating solid projected on a moving 
film. The general inverse problem is more complex and has yet to be solved for the general 
case. However, some conclusions taken from the direct problem can help the formulation 
of the inverse problem. An example is the case of regular polygons which was discussed 
in the previous section. 

Indeed, consider an image generated by the rotation of a regular polygon with con- 
stant angular velocity. Suppose also that the film moves with a constant speed. We want 
to know which polygon has generated that image. To solve this problem, suppose that the 
polygon rotates around its center of mass. If this centre is located at the origin of a system 
of Cartesian coordinates, then the shadow will be limited above and below by periodic 
curves. We can ensure, in advance, that if those curves are symmetrical with respect to 
the x axis then the polygon has a number of even sizes. If the curves are shifted, the 
polygon has an odd number of sides. To know the exact number of sides of the polygon it 
is enough to know the relationship between the values to, meaning the minimum distance 
from the superior curve to the x axis, and M, meaning the maximum distance from the 
superior curve to the x axis. The value of to is given by polygon's apotema and M is the 
distance of the center to the vertex of the polygon (see figure 19). 
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Figure 19: Image of a regular polygon 

7r / in \ 

Clearlyt 6 = — = arccos I — ) , where n is the number of edges of the regular polygon. 
n \M J 



Thus 



n 



arccos 



( m i 



Por instance, if we do not know in advance which the regular polygon generated the 
figure 15, by inspection of the image we may conclude immediately the regular polygon 
has a number of even sizes due the symmetry ot the upper and lower bounds. To find 
out the exact number of edges it is enough to calculate m and M. In fact, m = — and 

M = — — . Therefore, 



7T 



7T 



n 



arccos 



( a / 2 ) tt/4 
UV2/2J 1 



that is, the polygon is a square. 

For the case of figure 19, m - 
polygon is a equilateral triangle. 



1^3 



and M 



i-s/3 



. We conclude that n = 3 and the 



Note that the smaller the difference (M — m) , the greater the number of sides of the 
regular polygon. If m — > M then 



lim 



TV 



arccos 



( m ] 



+ OO, 



that is, in the limit, we have a rotating circle. 
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8. FINAL REMARKS. We have discussed in this short paper a simple problem. 
We believe however that it clarifies and illustrates some basic questions arising in direct 
and inverse problems. It becomes very clear, for instance, why convexity of domains of 
definition are so important for the uniqueness of the solutions, particularly for the case 
of inverse problems. 

The kinematic image can also be used as a tool to stimulate the abstract represen- 
tation of the shape of two-dimensional objects. In that sense it could be an interesting 
tool for educational purposes. Also more complex motion of the film and of the rotating 
object can be introduced. It is possible that the identification of intricate contours can 
take advantage of the comparison between two or more kinematic images corresponding 
to different motions. We believe that this problem is a starting point for several other 
setups including three-dimensional objects. 
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